Negative thermal expansion in ZnF 2 
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We have investigated temperature dependence of the lattice parameters and the unit cell vol- 
ume of ZnF2 by neutron diffraction and have discovered negative thermal expansion (NTE) at low 
temperature. To understand why this simple compound exhibits NTE we performed first principle 
calculations. These calculations reproduce qualitatively the experimental temperature dependence 
of volume. 

PACS numbers: 61.05.fm, 65.40.De 



The negative thermal expansion (NTE) in solids has 
attracted the renewed attention of condensed matter sci- 
entists ever since Sleight and coworkers^' discovered that 
ZrW 2 Os contracts over a wide temperature range of more 
than 1000 K. There are excellent review articles^bn the 
NTE of this type of so-called framework materials. How- 
ever NTE is known and has been studied experimentally 
and theoretically for a long time. Among these Si and Ge 
and other tetrahedrally bonded crystals at low tempera- 
ture are classic example^^. The NTE is however limited 
in the low temperature range in these materials, whereas 
in higher temperature range they exhibit normal positive 
thermal expansion. There exists a more general review 
article^ covering all types of materials that exhibit NTE. 
Here we report observation of NTE in diamagnetic ZnF2 
with the simple rutile-type structure. We also report the 
results of our ab-initio calculations that reproduce qual- 
itatively the observed NTE in ZnF 2 . 

The transition-metal difluorides MF 2 (M = V, Cr, 
Mn, Fe, Ni, Cu, Zn) with the rutile-type or distorted 




FIG. 1: (Color online) (Left panel) Rutile-type crystal struc- 
ture of ZnF2. The grey spheres are Zn and the green spheres 
represent F atoms. (Right panel) The stacking of ZnF6 octa- 
hedra in ZnF2. 
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FIG. 2: (Color online) Rietveld refinement of the ZnF2 crystal 
structure at T = 1.8 K. 



rutile-type crystal structure form an important class of 
materials with interesting magnetic and magneto-optic 
properties. In order to entangle magnetic effects from 
the lattice effects, the last of this series viz. the non- 
magnetic ZnF 2 has often been used to study the back- 
ground lattice effects. However even this non-magnetic 
solid showed anomalies in elastic constants at low tem- 
peratures not expected for ZnF 2 with no temperature- 
induced phase transition. The observed softening of 
C44 and C s = (Cu — Ci 2 )/2 at low temperature in 
ZnF 2 ha d bee n interpreted as an incipi ent fe rroelectric 
transitiorPlSI but has also been conteste d 11 * 12 !. 

Figure [I] shows schematically the rutile-type crystal 
strucure of ZnF 2 that crystallizes with the space group 
or P4/mnm. The unit cell is tetragonal with lattice 
parameter a = 4.7034, c = 3.1335 A at T = 296 K and it 
contains two formula units Z = 2. The two Zn 2+ ions are 
located at positions (0,0,0) and (1/2,1/2,1/2) whereas 
four F~ ions are located at (x, x, 0), (1— x, 1—x, 0), (1/2— 
x, 1/2 +x, 1/2), (1/2 + x, 1/2 -x, 1/2) with the positional 
parameter x — 0.303. The Zn 2+ ions are surrounded by 
six F~ ions to form slightly distorted octahedra. The 
octahedra are edge linked along the c-axis and corner- 
linked along < 110 > crystallographic directions. 

Neutron diffraction experiments were done on ZnF 2 
on the high intensity powder diffractometer D20 of the 
Institute Laue-Langevin in Grenoble. The 115 reflection 
from a Ge monochromator at a high take-off angle of 118° 
gave neutron wavelength of 1.868 A. Approximately 5 g 



2 




FIG. 3: (Color online) Temperature variation of the lattice 
parameters a, c, and the unit cell volume V of ZnF2 plotted on 
the left panels. The red curves in these figures represent the 
lattice parameter and the unit cell volume obtained by fitting 
the data by fifth order polynomials. On the right panels only 
the low temperature data are shown. 



ZnF2 powder samples was placed inside an 8 mm diame- 
ter vanadium can, which was fixed to the sample stick of 
a standard 4 He cryostat. We have measured the diffrac- 
tion intensities from ZnF2 as a temperature in the range 
1.8 - 320 K. The Rietveld refinementP of the diffrac- 
tion data was done by the Fullprof program^. The re- 
finement results from ZnF2 at T = 1.8 K is shown in 
Fig. [2] The agreement factors R (not corrected for back- 
ground) for pattern of this refinement were Rp — 3.13 



and R wp — 4.13. The corresponding conventional Ri- 
etveld R-factors were R p — 14.0 and R wp = 10.5. The 
goodness of the fit as given by \ 2 was X 2 = 2-07. 

Fig. [3] shows the temperature variation of the lattice 
parameters a, c, and the unit cell volume V of ZnF2 
on the left panels. The red curves in these figures rep- 
resent the lattice parameters and the unit cell volume 
obtained by fitting the data by fifth degree polynomials. 
Attempts to fit the data by Debye or Einstein functions 
in Griineisen approximation failed for a and V because of 
the negative thermal expansion at low temperaure. Fifth 
degree polynomial function fit the low temperature suc- 



FIG. 4: (Color online) (Upper panel) Temperature variation 
of AV/Vo- The blue circles are the data points and the dot- 
ted red curve is the DFT calculation of A V/Vo- (Lower panel) 
Temperature variation of the thermal expansion coefficient 
av = ^7 ^ • The blue and the dotted red curves correspond 
to the observed and calculated volume thermal expansion co- 
efficient av- The experimental values were determined from 
the differentiation of the fitted data with a fifth-order poly- 
nomial function. 



cessfuly. On the right panels only the low temperature 
data are shown. The lattice parameter a and the unit 
cell volume V exhibit minima at about 75 K. 

The upper panel of Fig. [4] shows the experimental 
normalized volume change given by 



AV _ V(T) - V 
~Vo : : V 



(1) 



where V(T) is the unit cell volume at temperature T and 
Vq is the volume at T = along with the values calculated 
within the density functional theory framework (DFT). 
The lower panel of Fig. [4] shows the experimental and 
calculated volume thermal expansion coefficients ay in 
the low temperature range given by 



a v (T) 



1 dV 
VdT' 



(2) 
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FIG. 5: (Color online) Temperature variation of the positional 
parameter x of the F atom and the two Zn-F bond distances 
dl and d2 of ZnF 2 . 



In order to check whether there exists any indication 
of incipient ferroelectric phase transitiorPQ2] in ZnF 2 at 
low temperature we refined the neutron powder diffrac- 
tion data by the Rietveld method and determined the 
positional parameter x and the two Zn-F bond distances 
di and g? 2 as a function of temperature. Fig. [5] shows 
these quantities. The absence of any anomalies suggest 
that apart from the negative thermal expansion (NTE) 



no further structural changes take place in ZnF 2 low tem- 
perature. 

We have done calculations to check whether we can re- 
produce NTE in ZnF 2 using first-principles DFT. Since 
DFT is a T = K approach, the finite temperature de- 
pendence is introduced within the framework of phonons 
by applying the so-called quasiharmonic approximation. 
The anharmonic effects are included uniquely via the vol- 
ume dependence of the phonon frequencies. For a set of 
volumes around the equilibrium one, the procedure con- 
sists of evaluating the total Hclmhotz free energy given 
by 



F(V,T) = E el (V) + F ph (V,T) 



(3) 



where E e j(V) and F p h(V,T) are the total ground-state 
temperature-free energy at constant volume as obtained 
directly from DFT and the phonon free energy extracted 
from subsequent lattice dynamical calculations for each 
volume, respectively. The free energy F(V,T) can be 
then used to evaluate thermodynamics of the material 
under study. 

The starting geometry for the calculations was the ex- 
perimentally refined ZnF 2 structure^. Relaxed geome- 
tries, total energies, phonon frequencies and volume- 
dependent phonon free energies were obtained using sim- 
ilar computational procedure described previously-ISHIll 

The underlying mechanism of NTE in ZnF 2 seems 
to be very similar to that for tetrahedral semiconduc- 
tors like Si, Ge, ZnS etc. with diamond and zincblende 
structures'"^. It is the excitations at low temperatures 
of the low-energy phonon modes with negative Griineisen 
parameters that are responsible for NTE in ZnF 2 . These 
modes are likely connected with the rigid-mode vibra- 
tions of ZnF 6 octahedra and their linkage along the a 
axis shown in the right panel of Fig. [T] 

In conclusion we have done neutron diffraction study 
of the temperature dependence of the crystal structure of 
the simple non-magnetic or diamagnetic transition metal 
difluoride ZnF 2 with rutile structure and have discovered 
negative thermal expansion (NTE) at low temperature. 
Our first principle calculations reproduces qualitatively 
this experimental result. 
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